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We present results of a lattice computation of the vector and axial- 
vector current matrix elements relevant for the semileptonic decay B° — ► 
p + l~&l. The computations are performed in the quenched approximation 
of lattice QCD on a 24 3 x 48 lattice at (5 = 6.2, using an 0(a) improved 
fermionic action. Our principal result is for the differential decay rate, 
dT/dq 2 , for the decay B° — > p + l~9[ in a region beyond the charm thresh- 
old, allowing a model-independent extraction of \V u b\ from experimental 
measurements. Heavy quark symmetry relations between radiative and 
semileptonic decays of B mesons into light vector mesons are also dis- 
cussed. 
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1 Introduction 



The experimental determination of the Cabibbo-Kobayashi-Maskawa (CKM) [[T], 
0| matrix elements is vital because they control the hadronic sector of the Stan- 
dard Model, determining CP violation and flavour mixing. The magnitudes and 
phases of these matrix elements must be established to make the Standard Model 
predictive, to determine the angles and area of the unitarity triangle, and to look 
for any inconsistency between the Standard Model and the experimental data 
which would point to new physics. 

The CKM element V u b is one of the most poorly known, currently uncertain to 
within a factor of two or three (at 90% CL) in magnitude ||. The determination 
of \V u b\ has traditionally been made from inclusive b — > u semileptonic decays ||, 
looking at the lepton energy spectrum beyond the endpoint for decays to charmed 
final states. These determinations rely on models incorporating nonperturbative 
QCD effects to relate the measured spectrum to theoretical predictions. Models 
used include the quark model || and bound-state models @-|| together with 
attempts to combine features from both ||. This results in the inclusion of an 
additional error reflecting the range of models used [ffO . 



In this paper we propose a model-independent method to determine |V^| 
from the exclusive semileptonic charmless B decay B° — > p + l~i?i, which should 
be measured with improved accuracy in B factories and at e + e~ and hadron 
colliders in the near future. Currently, only an upper bound exists [II] for the 
total decay rate, but new experimental results should be available soon f!2]| . 

We will show how lattice QCD calculations, which incorporate nonperturba- 
tive QCD effects in a systematic way, can be used to extract \V U b\ from experimen- 
tal measurements of B° — > p + l~i?i. Exclusive B decays have already proved useful 
for extracting the value of [V^l by studying the process B — > D*W[ []T3fl-[|r9||. 

We have also analysed relations between the exclusive processes B — > K*^ 
and B° — > p + l~vi which follow from heavy quark symmetry (HQS) P0| . We 
discuss the utility of these relations for determining \V u b\ from the experimental 
measurements of B(B —>■ K*j) and B(6 — > S7) f2"I||-|2"3 . 



2 Form Factors 

The matrix elements we will be considering are of the V — A weak current be- 
tween B and p mesons and of the magnetic moment operator between B and K* 
mesons p5 |. For B° — > p + l~vi the matrix element is, 



(p(k,rj)\u llM (l- l5 )b\B(p)) = v*%p, (1) 
with form factor decomposition, 

2V(q 2 ^ 



tub + m 



-e^spP^k - i{m B + m p )Ai(q )g fl/3 
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A 2 (q 2 ) . A(q 2 ) , . . 

+ * : (P + kj^qp-i — 2 2m p g A1 (p + A;)/?, (2) 

mg + m p g^ 

where g = p — k is the four-momentum transfer and 77 is the p polarisation vector. 
The form factor A can be written as 

A(q 2 ) = A (q 2 ) - A 3 (q 2 ), (3) 

where, 

2 m B + m p 2 m B - m„ , , , 

with Aq(0) = ^3(0). In the limit of zero lepton masses, the term proportional to 
A in equation @ does not contribute to the total amplitude and hence to the 
decay rates. Pole dominance models suggest that V, A4 for % = 1,2,3 and A 
correspond to 1~, 1 + and 0~ exchanges respectively in the t-channel @. 

The main contribution to the B — > K*^ decay comes from the matrix element 

(K*(k, V )\sa^q»b R \B(p)) = £c;T 4 (g 2 ), (5) 

i=i 

where q = p — k as above, rj is now the K* polarisation vector and 

Cj = 2e^ XpV »p x ko, (6) 

cl = v^ 2 B-^K*)-v-q{p + k)^ (7) 

CJ = ^-gL- g ' (p + *)„)■ (8) 

For an on-shell photon with q 2 = 0, T 3 does not contribute to the B — > if* 7 
amplitude and 7\ and T 2 are related by, 

T 1 (g 2 =0)=zT 2 (g 2 =0). (9) 

Hence, for B —>■ if*7, we need to determine T\ and/or T 2 at the on-shell point. 

Neglecting corrections suppressed by inverse powers of the heavy quark mass 
M, the following relations hold when g 2 is close to the maximum recoil value 

ql sx = {M-m p , K *Y 



VQ/y/M = const, A^y/M = const, A 2 Q/y/M = const, (10) 

where arises from the leading logarithmic corrections and is chosen to be 1 at 
the B mass IE? 
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In the calculations reported below, we will use (5q = 11 in the quenched approx- 
imation and Aqcd — 200 MeV. The matrix elements of the V — A current and 
the magnetic moment operator between B and identical^ light final-state vector 
mesons of mass m are related in the infinite heavy quark mass limit according to: 



V(q 2 ) = 2T x (q 2 



A l {q 2 ) = 2iT 2 {q 2 



(12) 



for values of q 2 not too far from g max , or equivalently, for uj close to 1, where 



CO 



v ■ v 



M 2 + m 2 - q 2 
2Mm 



(13) 



Here, v and v ' are the four-velocities of the heavy-light pseudoscalar meson and 
the light vector meson respectively. 

The equations in (|TJ) relate the physical decay processes B° — > p + l~vi and 
B~ — > p~j. If light flavour 577(3) symmetry is respected then equation (|i~2j) 
relates the processes B° — > p + l~vi and B — > K*^ ||2U|| . In the lattice calculations 
reported below we will test equation (|12"D directly for different values of the heavy 
quark mass using identical light meson states for both matrix elements. 

In the literature, versions of equation ( |T2| ) appear which rely on lowest order 
HQS but incorporate all 1/M corrections from kinematics p0| , [21], ^7]0. The 
modified relations become: 



m (q 2 



2iT 2 (q 2 



q 2 + M 2 -m 2 V{q 2 ) M + m 



2M M + m + 2M 



(M+m) 2 - q 2 (M — m) 



-Mq 2 



V(q 2 



(14) 



2M(M 2 



M + m 



+ 



q 2 + M 2 - m 2 A x { 



2M 



M -m 



(15) 



3 Lattice Details 

Lattice calculations with propagating quarks provide matrix elements for heavy 
quarks around the charm mass over a range of q 2 straddling q 2 = 0. In extracting 
form factors from these matrix elements, we can reach g^ ax for A\ and T 2 only, 
because the coefficients determining the contribution of the other form factors to 
the matrix elements of equations (0) and ([5]) vanish at this kinematic point. To 
obtain results relevant for B decays, we need to extrapolate in the heavy quark 
mass M to the b quark mass. This is simple for fixed uj, but, at the B scale, 
produces a range of q 2 values near g^ ax and far from q 2 = 0. 

3 We assume isospin symmetry between the light u and d quarks. 

4 To incorporate properly the dynamical 1/M corrections it is necessary to consider operators 
of dimension four in the heavy quark expansion in addition to those of dimension three occurring 
at leading order. 
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Figure 1 Labelling of quark hopping parameters for three-point correlator calculation. 



The results described below come from 60 £77(3) gauge configurations gener- 
ated by the UKQCD collaboration on a 24 3 x 48 lattice at (3 = 6.2 in the quenched 
approximation. The 0(a) improved Sheikholeslami-Wohlert (SW) |28| action 



was used for fermions, with "rotated" fermion fields appearing in all operators 
used for correlation function calculations p9| . The inverse lattice spacing deter- 
mined from the p mass is a -1 = 2.7(1) GeV [30]. Other physical quantities will 
lead to slightly different values for the lattice spacing (a" 1 = 2.7-3.0 GeV 
The scale uncertainty should be reflected in the results for dimensionful quanti- 
ties. 

Three-point correlators of the heavy-to-light two-fermion operators with a 
heavy pseudoscalar meson (the U B" meson) and a light vector meson were cal- 
culated, as illustrated in figure p]. Matrix elements were extracted from these 



correlators by the method detailed in [32]- 35 . Four heavy-quark hopping pa- 



rameters, Kh = 0.121,0.125,0.129,0.133, were used. For the propagator con- 
necting the current operator to the light meson operator, two kappa values, 
K a = 0.14144,0.14226, were available. The subscript a is for "active" to con- 
trast this propagator with the "spectator" propagator joining the heavy-light 
meson to the light meson. These n a values straddle that for the strange quark, 
0.1419(1) |30| . For n h = 0.121,0.129, we used three light spectator hopping pa- 
rameters, ki = 0.14144,0.14226,0.14262, and for n h = 0.125,0.133 we used Ki = 
0.14144 only. The critical hopping parameter at this (3 is K crit = 0.14315(1) 0. 

The lattice calculations were performed with the heavy meson spatial momen- 
tum of magnitude or 1, in lattice units of 7r/12a. The momentum injected at 
the operator insertion was varied to allow the modulus of the light meson spatial 
momentum to take values up to y/2 in lattice units (although some of the mo- 
mentum choices were too noisy to be used in fits). We refer to each combination 
of heavy and light meson three- momenta as a channel with the notation |p| — > |k| 
in lattice units (for example 0— > 1 or 1 - » lj_ where the subscript _L indicates 
that p and k are perpendicular). 

The results below have been obtained using uncorrelated fits for the extrap- 
olations in the heavy quark mass. The extraction of the form factors from the 
three-point correlation function data used correlated fits fl32| , |33fl . Statistical er- 
rors are 68% confidence limits obtained from 250 bootstrap samples. 
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To make the best use of HQS, we pick momentum combinations which keep 
lo constant as the heavy-light meson mass varies. This allows us to scale in 1/M 
from the charm to the bottom mass scale and gives the form factors for the B 
decays as functions of u, and therefore as functions of q 2 . When the heavy-light 
meson is at rest, uj is independent of the heavy-light meson mass. As shown in 



reference |36[ there are some additional channels, where the heavy-light meson is 
not at rest, but p and k are perpendicular, for which uj is very nearly constant 
as the heavy meson mass varies, and which can be used for heavy quark mass 
extrapolations. For those channels where uj is not strictly constant but which we 
nevertheless use in our study, the difference between the average uj value used in 
our analysis and that for each value of the heavy quark mass is always less than 

3% g. 

For relating lattice results to the continuum, we have used the perturbative 
values for the renormalisation constants of the vector, axial and magnetic moment 



operators pTfl : 



Z v = 0.83, Z A = 0.97, Z a = 0.98. (16) 

The SW improved action reduces the leading discretisation errors from 0(a) 
in the Wilson fermion action to 0(a s a), but for quark masses tuq around that 
of the charm quark, a s rriQa can be of order 10%. An estimate of the lattice arte- 
facts can be obtained by comparing values of lattice renormalisation constants, 



computed from different matrix elements, in a non-perturbative way 38 . Nu- 



merical calculations have confirmed that errors of order 10% are present at our 



value of f3 in the matrix elements of vector and axial vector currents [15, 16]. We 
will therefore allow for an extra 10% systematic uncertainty in our results due to 
possible discretisation errors. 



4 Semileptonic and Radiative Decays of Heavy- 
Light Mesons and Heavy Quark Symmetry 

In this section we study the relations between the semileptonic and radiative 
decay processes of heavy-light mesons to light vector mesons. We will check the 
HQS relations of equation (|12|) and the size of 1/M corrections to these relations 
for different values of the heavy-light meson mass. For testing these relations we 
use our most accurate light-quark data, Ki = n a = 0.14144: this ensures that the 
light degrees of freedom are the same for the semileptonic and radiative processes 
and allows us to check directly, without corrections due to SU(3) symmetry 
breaking, the validity of HQS as M increases. 

Specifically, we compare the form factors A\ and V for the pseudoscalar to 
vector semileptonic decay, 

P hl { Kh ; ki = 0.14144) -> VuiKi = 0.14144)/^, (17) 
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with the form factors T\ and T 2 for the pseudoscalar to vector radiative decay, 



Phii^h] Ki = 0.14144) -> V n (Ki = 0.14144)7, (18) 

for different values of the heavy quark mass, Kh = 0.121, 0.125, 0.129 and 0.133, 
and for different values of to close to u — 1. The corresponding heavy-light pseu- 
doscalar masses are in the range 1.6 GeV to 2.5 GeV. We have also extrapolated 
to the B scale and to the infinite heavy quark mass limit. We have used five 
momentum channels in our analysis^ -»• 0, -»• 1, -»• \/2, 1 -»• and 
1 -> U- 

In all cases the ratios V/2T\ and Ai/2iT 2 are consistent with unity in the 
heavy quark mass limit, as predicted by HQS, although for certain channels the 
ratio V/2Ti departs from 1 by up to 75% at the charm scale. This constitutes 
a non-trivial test of HQS predictions. These results are illustrated in figures [| 
and p. 

In figure |^ we show the ratios V/2Ti (for momentum channels — >• 1 and 
1 — >• 0) and Ai/2iT 2 (for momentum channels 0^0 and 1 — > for different 
pseudoscalar meson masses, allowing for linear and quadratic 1/M corrections 
to the infinite mass limit predictions of equation (0). The linear and quadratic 
extrapolations in 1/M agree within errors for the extrapolated values, with the 
exception of three cases of the infinite mass limit points (one of which is shown 
in figure 0). We trust the linear extrapolations more than the quadratic ones: 
we are fitting to four points only in a limited region of inverse heavy-light pseu- 
doscalar mass, and quadratic fits can amplify accidental quadratic effects in the 
fitted points. All results in the remainder of this section will refer to linear 
extrapolations. 

Figure |3] shows the ratios Vj2T\ and Ai/2iT 2 for five values of u at three 
different heavy-light pseudoscalar masses, around the D mass, around the B 
mass and in the infinite mass limit. The results around the D mass are our 
measured values at = 0.129 (which corresponds very closely to the charm 
quark). The B mass and infinite mass limit results are extrapolations. The HQS 
predictions of equation (|12"D are well satisfied for both ratios in the infinite mass 
limit. The ratio V/2T\ shows large 1/M corrections of the order of 75% at the 
D scale and 20% at the B scale. Corrections of about 30% at the D scale and 
about 10% at the B scale were previously observed for the pseudoscalar meson 
decay constant fp J^J and for the form factor hy in heavy-to- heavy 0~ — > 1~ 
semileptonic decays |16|Q- m contrast, the ratio Ai/2iT 2 exhibits small 1/M 



corrections even at the D meson scale. A similar situation occurs for the form 
factors h + and h& x in heavy-to-heavy 0~ — > 0~, 1~ semileptonic decays fL5" , p!6[f|. 



5 ^4i(^ ax ) and T2(q^ nax ) can be determined directly from the — > data, but V(q'^ nax ) and 
^i(<Zmax) have to be obtained by an extrapolation in q 2 from the measured data. 

6 Short distance corrections have been accounted for in the form factors hv, h+ and Kai- 
7 Note that h + and liA-t are protected from 1/M corrections at zero recoil by Luke's theo- 
rem and the leading corrections are of order 1/M 2 . 
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Figure 2 V/2T\ and A1/21T2 as a function of the ratio of the light vector meson 
mass over the heavy-light pseudoscalar meson mass, for different heavy quark masses 
(squares) and various momentum channels. For the 0^0 channel u> is 1 exactly. For 
the remaining channels, the value of uj shown in each plot has an error of less than 1 
in the last digit. We show both linear (solid) and quadratic (dashed) extrapolations in 
1/M. The extrapolations to the B scale and the infinite mass limit are indicated by 
diamonds and crosses respectively. 



Inclusion of kinematic 1/M factors, as shown in equations ( l4|) and (la|) re- 
duces the size of the 1/M corrections by about a factor of two for the relation 
between T\ and V, giving results in agreement with those in reference but 
has little effect on the T2 and A\ relation (note that the Ti~A\ relation in equa- 
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Figure 3 Ratios V/2T\ and A\/2iT2 for five values of oj at three different heavy-light 
pseudoscalar masses, around the D mass (crosses), around the B mass (diamonds) 
and in the infinite mass limit (squares). The horizontal solid line denotes the HQS 
prediction in the infinite mass limit. 



tion ( P~5| ) reduces to the result in equation ( |T2D at to 



5 Cabibbo Suppressed Decays 

For this exploratory study we did not have a complete set of light quark kappa 
values for every heavy-quark kappa value, as discussed in section ||. This pre- 
vented us, in [36|, from performing a reliable chiral extrapolation to obtain results 
for the semileptonic decay B° — > 7r + /~z/;. Here we are concerned with semilep- 
tonic decays with a light vector meson in the final state: in contrast to the case 
of a pion, which is a pseudo-Goldstone boson, we do not expect a sizable effect 
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Figure 4 A\/Za for two different momentum channels and six different combinations 
of spectator and active quark masses as a function of the inverse spectator quark 
kappa value. Squares denote n a = 0.14144 and diamonds denote n a = 0.14226. The 
diamond points are slightly displaced horizontally for clarity. The heavy kappa values 
are Kh = 0.121 for the — * 1 channel and 0.129 for the 1 — » channel. The vertical 
dashed line marks the chiral limit for the light spectator quark. 



from the chiral extrapolation. 

Indeed, in figure ^ we show the dependence of the form factor A\ on the active 
and spectator light-quark masses, for two different combinations of momentum 
channel and heavy quark mass, = 0.121 and 0.129, for which six combina- 
tions of light quark kappas were computed. It can be seen that in both cases Ai 
remains practically constant as the light spectator quark mass decreases (hori- 
zontal movement on the plots) but there is a dependence on the active quark 
mass (vertical movement). Note that the value of u, and hence q 2 , depends on K[ 
as well as K a: and therefore the results of figure [| indicate that Ai' ,Ka (oj(ku K a)) 
is nearly independent of Ki, in which case, 

A^ M (u(K cli t,K ait )) « Ar ai4144 ' Kcrit (w(k,=0.14144, « crit )). (19) 

Similar independence of the spectator mass was found in the study of the form 
factor T2 in the radiative decay B — ► K*j . As we have shown in the previous 
section, A\ is equal to UTi to good approximation for different heavy quark 
masses and different momentum channels and therefore the results of give us 
further evidence of the spectator mass independence of A\. 

The situation is less clear for the form factors A2 and V owing to the larger 
statistical errors. A 2 and V appear to follow the same pattern as A±, but we 
cannot dismiss a possible mild dependence on the light spectator mass. In figure [5] 
we show two examples of the dependence of the form factors V and A 2 on the 
active and spectator light quark masses. 
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Figure 5 F(g^_ >0 )/Zy and ^(gfj— >i)/^4 f° r ^ wo different momentum channels and six 
different combinations of spectator and active quark masses as a function of the inverse 
spectator quark kappa value. Squares denote n a = 0.14144 and diamonds denote 
K a = 0.14226. The diamond points are slightly displaced horizontally for clarity. The 
heavy kappa value is Kh = 0.129 in both cases. The vertical dashed line marks the 
chiral limit for the light spectator quark. 



We are going to assume the form factors are independent of the light spectator 
mass, and therefore we will extrapolate to the chiral limit only for the active 
light quark, as was suggested for A\ in equation fll9|) . This will allow us to 
use results for all four heavy quark masses to guide our extrapolations to the B 
mass. To quantify the systematic error induced we have compared, for the two 
heavy kappa values for which we have six combinations of light quark masses, 
the chiral-extrapolated form factor values for each momentum channel using two 
procedures: 

1. Using all six data points we can extrapolate to the chiral limit for both the 



spectator and active quarks (see reference |33 



2. Fixing the spectator kappa value to 0.14144 and assuming complete inde- 
pendence of the form factors on the spectator quark mass, we perform the 
chiral extrapolation for the active quark only, as indicated in equation (|19|) 
for A 1 . 

For A\ the two procedures give the same results within 5% except for those chan- 
nels involving a momentum of \[2 in lattice units for the light vector meson, 
where the variation can be 10%. These errors are comparable with the statistical 
ones. For procedure 1 above, the chiral extrapolation for the \/2 channels is less 
reliable than for channels with lower light meson momentum, because the sta- 
tistical fluctuations for three-point functions with the lightest light-quark masses 
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channel q 2 / GeV 2 A l A 2 V 



- 


-> 


20.3 


0.461 2 . 






- 


-> 1 


17.5±1 


0.4312 


0.8^ 


1.611 


- 




15.311 


0.3911 


0.7l? 


1.211 


1 - 


-> 


19.711 


0.4611 






1 - 


+ i± 


16.7l| 


0.3811 


0.611 


1.511 


1 - 


-> v^± 


14.411 


0.391^ 


0.711 


1.412 



Table 1 Values of the form factors for B° — > p + l~i>i for different values of g 2 close to 
Qmax- For the 0^0 channel, w is 1 exactly, so that g 2 is fixed up to (tiny) experimental 
errors in the physical meson masses. For A 2 and V, the zero recoil channel, — > is not 
measured. In addition, the large statistical errors in A2 and V for the 1 — ► channel 
prevent a reliable extrapolation in 1/M. Quoted errors are purely statistical. Adding 
systematic errors from spectator quark flavour symmetry breaking and discretisation 
in quadrature produces a further 11% error in Ai, 20% in A2 and 15% in V. 



(ki = 0.14226,0.14262) are larger. These fluctuations can account for an impor- 
tant part of the 10% variation observed between the two procedures. We have 
decided to use 5% to estimate, for all momentum channels, the systematic error 
on Ai induced by this assumption of independence of the light spectator quark 
mass. 

The statistical fluctuations are much larger for A 2 and V, making it diffi- 
cult to check the difference between the results from the two chiral extrapolation 
procedures above. We will adopt a conservative position and will admit a 20% 
systematic error in A 2 and a 10% error in V, based on the maximum discrep- 
ancies observed between the two procedures. These discrepancies are within the 
statistical error bars. 

In figure || we show extrapolations to the B mass scale for the form factors 
Ai, A 2 and V, allowing for linear and quadratic 1/M corrections to the infinite 
mass limit predictions of equation (f[0|). The extrapolated values agree well for 
the linear and quadratic fits: we will quote results only from the linear fits in 
the following, unless stated otherwise. We have used six momentum channels: 
the same five channels used in the previous section together with the channel 
1 — > y/2j_. For A 2 and V the momentum channels — * and 1 — >• are 
extremely noisy and we cannot control the extrapolations!. The measured values 
of the form factors, extrapolated to the B scale for different values of q 2 close to 
g^ax, are listed in table |TJ. 

In figure 0a we show Ai(q 2 ) at the B scale, together with three fits to its q 2 

8 The — > momentum channel is not measured directly for A2 and V . 
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Figure 6 Linear (solid line) and quadratic (dashed line) 1/M extrapolations of A±, A2 
and V in selected momentum channels as functions of m p /M. Squares denote measured 
points. Diamonds and crosses mark the extrapolations to the B scale for the linear 
and quadratic fits respectively. 



dependence: 



r^i(o) 

A 1 (0) 



AM 2 ) 



<l 2 / m lo\e 

Mo) 



y (1 - q 2 /m 2 dipolc ) 2 



constant, 
pole, 

dipole. 



(20) 



Table |2| shows the fit parameters for each q 2 dependence, with the corresponding 
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Figure 7 (a): form factor Ai(q 2 ) for the decay B° — » Squares are measured 

data extrapolated to the B scale at fixed uj. The three curves and points at q 2 = 
have been obtained by fitting the squares using the q 2 dependences from equation (pc|): 
constant (dashed line and octagon), pole (solid line and diamond) and dipole (dotted 
line and cross). The point at q 2 = for the dipole fit has been displaced slightly for 
clarity, (b): linear (solid) and quadratic (dashed) extrapolations of Ai(0) in 1/M to 
the B scale according to equation (|2l]) . Squares are obtained from pole fits to the 
measured data for the different heavy quark kappa values. The diamond and cross give 
the form factor at the B scale. 



fit type 


am 


"Vie/ GeV 


"idipoio/ GeV 


X 2 /dof 


constant 

pole 

dipole 


0.43+1 
0.2711 
0.25±| 




9±j 


1.3 
0.2 
0.2 



Table 2 Fit parameters for different q 2 dependences for the form factor A\(q 2 ) at the 
B scale. Quoted errors are statistical only. 



X 2 per degree of freedom. 

The figure and table indicate that constant-in-g 2 behaviour for Ai does not 
fit the data well. This feature is even more evident in lattice studies at the D 
scale, where fits to constant behaviour are poor |33|-|35]], (42], [43]]. Our data 
favour a pole type behaviour with a pole mass in the expected range for a 1 + 
bu resonance ^3|. Dipole (and in general higher powers: tripole, . . . ) and pole 
fits are hardly distinguishable in the physical range < q 2 < q^^. In a dipole 
fit, the mass parameter mdi po ie is roughly given by mdi po i e ~ V^m P oic and hence 
pole and dipole fits agree in the limited range of values of q 2 , well below ^p 0le , 
explored in figure 0a. An alternative procedure is to use pole fits for extracting 
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Ai(Kh, q 2 =0) from our measured data, at different heavy quark masses around 
that of the charm and then to extrapolate Ai(0) in 1/M to the B scale, as shown 
in figure ^b, by using |32], |36], pBf 

Ai(0)6M 3/2 = const(l + j/M + 5/M 2 + ■■■). (21) 

This gives a result for Ai(0) at the B scale (using linear or quadratic extrapola- 
tions in 1/M) of, 

, . 2 n . (0.18 ±0.02 linear , no . 
A l( g 2 =0;m B ) = ( a22 ^ quadratic - (22) 

The quoted errors are purely statistical. A further 11% systematic error should 
be added, obtained by combining in quadrature the systematic errors from the 
spectator quark flavour symmetry breaking and discretisation effects. Incorporat- 
ing the systematic error makes this method of extracting A^O) consistent within 
errors with the method presented in table |2|. 

In conclusion, we find that A\(q 2 ) at the B scale is fitted by a single pole 
form, with the parameters given in table 0. The errors in table |2| are statistical 
only: as mentioned above, a further 11% systematic error should be added. In 
table fj we compare our results with other theoretical calculations: the agreement 
is generally good. Previous lattice results W\ rely on the assumption of pole 



behaviour using a lattice determination of the bu 1 + resonance (B\) mass and 
the value of Ai(q 2 ) at a single q 2 value. In contrast, in this paper, we have tried 
to determine the q 2 dependence of A\. 

Statistical and systematic uncertainties in the present study have prevented 
us from extracting the q 2 dependence of the A 2 and V form factors. 

Much effort has recently been devoted by the lattice community to deter- 
mining the q 2 behaviour of T\ and T 2 around the B scale, without producing a 
definitive conclusion. The q 2 dependence found here for A\ and the results from 
the previous section for the ratio A\jHT2 may support a pole- type behaviour 
for T 2 at least in a region around g^ax- Sum-rules calculations [[11], ^5], |46|] and 



a recent calculation by B Stech EDI find that V has a more pronounced q 2 de- 



pendence than Ai, which is consistent with having dipole type behaviour for V 
and pole behaviour for A%. In the previous section we found that around the 
B scale T\ is roughly equal to V. This supports a dipole behaviour for T\ and 



consequently further favours a pole behaviour for T 2 [32, 36, 14 . 



6 Extraction of |V^| 

In this section we employ our lattice determination of the form factors to calcu- 
late the differential decay rate dT/dq 2 for the decay B° — ► p + l~vi in the region 
near g^ ax . Experimental data for this region will enable a model- independent 
determination of \V u h\. 
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Reference 



42 1 



This work 
ELC "a" 
ELC "b" 
APE "a" 
APE "b" 
Sum rules ET 



32 1 



43 1 



Sum rules [45 



Sum rules [46 



Quark model || 
Quark model 
Quark model fl47 | 
HQS & xPT |H 



^i(O) 
0.27+1+1 
0.25 ±0.06 
0.22 ±0.05 
0.29 ±0.16 
0.24 ±0.12 
0.24 ±0.04 

0.5 ±0.1 
0.35 ±0.16 
0.05 
0.28 
0.26 ±0.03 
0.21 



Table 3 Values for A^O) for B° -> /? + /~^ from various theoretical calculations. The 
second set of errors in our quoted value denotes our estimate of systematic effects from 
the spectator quark flavour symmetry breaking and discretisation. 



It has previously been suggested that determinations of the decay rate for 
B° — > p + l~vi near q 2 = 0, combined with the experimentally-measured branch- 
ing fractions B(B — ► K*j) and 5(6 — > s'j) would also allow a determination of 
\V u b\. However, one method relies on three different measurements, B(B — > K*^), 
B(b — > 57) and gT(7? — > p + l~vi)/dq 2 \ q 2 =0 , together with the theoretical determi- 



nation of the ratio T x ~* 7 (0)/7L ^ p (0) p2L 23]. The evaluation of this ratio of 



form factors from lattice calculations involves difficult and model-dependent ex- 
trapolations from the vicinity of g^ ax to q 2 = 0. Another method uses the ratio of 
T(B —>■ K*j) and lim^^o dT(B° — > p + l~i?i)/dE p dEi evaluated on a curve where 
q 2 = 4Ei(rriB — E p — E{) [^1|], which is independent of any hadronic form factor. 
However, it relies on the validity of the leading order HQS relations of equa- 
tions ( |I~4D and ( [IB] ) between the form factors for the radiative and semileptonic 
decays in the region around q 2 = 0, far from the zero recoil point. Furthermore, 
light flavour £77(3) symmetry is used to relate the K* and the p mesons. Correc- 
tions to these approximations will induce systematic errors in the determination 
of \V u b\. For example, the sum rule calculation of the relevant form factors in 
reference [51] finds about 7% corrections to the leading HQS relations and about 
25% corrections from SU (3) breaking. 

We propose to look at the region around g^ ax for the process B° — ► p + l~i> h 
beyond the region of charm production which complicates the experimental de- 
termination of b — > u transitions at low q 2 (note that \V u b/V c b\ 2 determined from 
inclusive decays is expected to be less than 0.01). The methods referred to above, 
based on determinations of the decay rate near q 2 = 0, suffer from this experimen- 
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tal difficulty in addition to possible theoretical uncertainties. Lattice techniques, 
in contrast, allow the form factors to be measured directly near q 2 ^^, using only 
a 1/M extrapolation motivated by HQS, and avoiding the problematic model- 
dependent extrapolation to q 2 = 0. This model-dependence currently plagues 
the lattice determination of B(B — > K* r f) which requires knowledge of the form 
factors Ti and/or T 2 at q 2 = 0, where they are equal. 

The differential decay rate for B° — ► p + l~vi is given by 



2 \\r . |2 



dT(B° -»■ p+l-ih) G 2 F % 



ub\ 2rw 2 



dg 2 1927r 3 m| 



q 2 [X(q 2 p 



H+(q 2 )\ 2 +\H~(q 2 )\ 2 + \H (q 2 )\ 2 ) (23) 



where A(g 2 ) = (m 2 B + m 2 — g 2 ) 2 — 4m 2 3 m 2 . if comes from the contribution of 
the longitudinally polarised p and is given by 

H\q 2 ) = tt^{« -m 2 p - q 2 ){m B + m p )A x tf) - ^^A 2 (q 2 )}, (24) 
2m p Vf I m B + m p J 

where k is the momentum of the p in the 5-meson rest frame. if ± correspond 
to the contribution of the transverse polarisations of the vector meson and are 
given by 

H±(q 2 ) = -Um B + m p )Aitf) =F 2m ^' k ' T/( g 2 )j, ( 25 ) 
I (m B + m p ) J 

Looking for semileptonic decays beyond the charm production threshold (de- 
termined by the semileptonic B — ► Z?ZPj decay) will make sense only if there 
are enough events. In figure || we show the differential decay rates dT/dq 2 for 
B° — > p + l~vi and 5° — ► 7r + / _ z/;. For 5° — > p + l~vi we assume that Ai, A2 and 
V are all given by single pole forms with a common pole mass of 5.3 GeV and a 
normalisation determined by our measurements for the — ► 1 momentum chan- 
nel. For B° — > 7T + l~i?i we use the results of reference |36[ with a dipole form for 



the form factor / + with / + (0) = 0.24 and a mass parameter mj+ = 5.7 GeV. 
These curves are meant to be illustrative, based on reasonable hypotheses: for 
B° — > 7T + l~i?i the results in [3B| were obtained for unphysically large u and d 



quark masses; for B° — > p + l~vi we do not necessarily believe that the three form 
factors involved are all simultaneously single pole types with the same pole mass 
(the pole mass of 5.3 GeV is an estimate of the average of the 1~ and 1 + bu 
resonance masses). The points we wish to make are: 

1. For B° — > p + l~vi the events near q 2 = have already been recommended for 
study, but the expected number of events at large q 2 is at least comparable. 

2. In the region beyond the charm threshold, the rates for B° — > p + l~vi and 
B° — > 7T + l~Pi are at least comparable (B° — ► 7r + /~z// is not ten times more 
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Figure 8 Differential decay rates for B° — > p + l^vi and B° — ► 7r + /~P;. The data 
points and solid curve are for B° — ► p + l~9i: squares have contributions from our 
measured values for A±, A2 and V, crosses have A\ only, and the crosses are offset 
slightly for clarity. The dashed curve is for B° — > 7r + /~z^. The curves are illustrative, 
as described in the text. The vertical dot-dashed line marks the upper limit in q 2 for 
charm production. 



common than B° — > p + l~i?i, for example), so both should be sizable con- 
tributions to the inclusive b — > u event rate and either or both can be used 
for the determination of \V u b\. 

The lattice calculations of the form factors in the previous section provide 
data points for dT/dq 2 over a sizable fraction of the region in q 2 above the charm 
endpoint, as can be seen in figure ||. The uncertainties in determining A 2 and V 
in the — > and 1 — > momentum channels do not prejudice our measurement of 
the differential decay rate. The rate is forced to vanish for kinematical reasons at 
<?max so t ne 0^0 channel form factor values of A2 and V are irrelevant. Moreover, 
the contribution to the differential decay rate of these two form factors in the 
1 — >■ channel is highly suppressed (the 1 — > channel has a q 2 near g^ ax and the 
A2 and V contributions are suppressed relative to that of A\ by a factor k 2 /m 2 3 as 



shown in equation (23)). Over the range of q 2 for which we have measurements, 
the differential decay rate is dominated by the form factor A\. This is shown by 
the difference between the points marked by squares and crosses in figure ||. 

In table |4] we give the numerical values we have obtained for the differential 
decay rate at several q 2 values. We have parametrised the results of the table by 
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q 2 


dT/dq 2 Jjr* 


dq 2 dT/dq 2 


20.3 
19.71} 
17.5±| 
16.711 
15.31^ 
14.41^ 


0.0 

0.19l| 
0.571^ 
0.581^ 
0.61} 
0.8l 2 


0.0 
0.081} 
0.91} 
1.31} 
2.3l| 
3.01^ 


charm threshold 
11.6 




5-411 



Table 4 Differential and partially integrated decay rate for B — * p + l v\ measured 
on the lattice for various q 2 values, q 2 is given in units of GeV 2 and T in units of 
|K&| 2 10 -12 GeV. For the ^ channel, u> is 1 exactly, so that q 2 is fixed up to (tiny) 
experimental errors in the physical meson masses. The integration has been performed 
using the parametrisation of equation (p6|) . Quoted errors are purely statistical — a 
further 24% systematic error should be added as discussed in the text. 



fitting our measured points for dT/dq 2 to the form 

W^ driS ° 7f'~ P,) = mh*MW + ^ - O), (-) 

obtained by retaining the phase space dependence and expanding H ± and H° 
from equation (|23|) around g^ ax . The results of the fit are 

a 2 = 2111 GeV 2 , 
b = (-8l^)10- 2 GeV- 2 , (27) 

where the quoted errors are statistical. A 24% systematic error should be added 
as discussed below. 

With this parametrisation, we can partially integrate the differential decay 
rate, giving the results shown in the last column of table [|. The table also includes 
the partially integrated decay rate from the charm threshold to q^ax, although 
higher order terms in the Taylor expansion in equation (|26|) could become sizable 
at the lower end of this range, where we do not currently have measured points. 
The results of table f| and equation (|27j ) have been obtained by extrapolating 
the form factors linearly in 1/M. Quadratic extrapolations in 1/M give results 
differing in the last significant figure, which is always well within the statistical 
errors. In figure |^ we show our measured data together with the fit, including 
68% confidence level bounds. 

Comparison of any of the values in table [| with experimental measurements 
will allow an extraction of \V u b\ with less than 10% statistical and 12% systematic 



19 




Figure 9 Differential decay rate as a function of q 2 for the semileptonic decay 
B° —* p + l~vi. Squares are measured lattice data, solid curve is fit from equation (p6|) 
with the parameters given in equation Dashed lines are 68% confidence limits 

(statistical errors only). The vertical dotted line marks the charm threshold. 



theoretical uncertainties. Another way to determine \V u b\ uses the parametrisa- 
tion of equation (^) and our results of equation (0). It parallels the method used 
to extract \V c b\ from semileptonic B — > D* decays (see references [[L3|l -||19||). The 
factor a 2 (l + b(q 2 — q^ax)) i n equation (|26|) parametrises long distance hadronic 

dynamics (it is the analogue of the function £, 2 (uj) in for the B — > D* case). 
The idea is to use experimental measurements of dT/dq 2 extrapolated to g^ ax to 
extract a value for iK^^a 2 and compare with our theoretical determination of a 2 . 
Here a plays the role of £(1) in the B — > D* extraction of \V c b\- 

The lattice simulation provides a systematic determination of the overall nor- 
malisation of hadronic dynamical effects, parametrised by a. For the extraction 
of \V c b\ from B —>■ D* decays HQS provides such a normalisation (the Isgur- 
Wise function is 1 at zero recoil, up to 1/M 2 and short distance perturbative 
corrections). For the heavy-to-light B° — > p + l~vi semileptonic decay, such a 
normalisation is not provided by HQS. 

In the fit to equation (p6|) the — > channel is not used because dT / dq 2 \ q 2^ 
is zero and our fitting function automatically incorporates this feature. Our value 
for ^4i(q , max) fr° m the — > channel therefore gives an independent measurement 
of a 2 because V and A 2 do not contribute at this kinematic point. We find, using 
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the value from table ^ and quoting statistical errors only, 



a 2 = S(m B + m p YAi(q 2 max ) = 23t* GeV\ (28) 

in excellent agreement with a 2 determined by the fit to the other channels. This 
suggests that discretisation and other systematic errors have not conspired to 
change drastically the shape of the overall q 2 dependence. Our final result will 
use the value of a 2 , given in equation (p7|) , determined from the fit to these other 
channels, because this does not rely on a single measurement. 

To estimate systematic uncertainties, we have considered four possible sources 
of errors: quenching, determination of the lattice spacing in physical units, light 
spectator mass independence of the form factors and discretisation errors. 

The exact effect of ignoring internal quark loops is an unknown systematic 
error. However, quenched lattice calculations of form factors for D — > K, K* 
semileptonic decays (see and references therein) give results in agreement with 



world average experimental values, while quenched calculations of B — > D, D* [ITS 



decays have allowed extractions of \V c b\ in agreement with other determinations 
of that quantity. This gives us some confidence that errors due to quenching are 
most likely within the statistical errors for processes involving a heavy quark. 

The error arising from uncertainties in the value of the lattice spacing in 
physical units should be minimised in our calculation, because we have evaluated 
only dimensionless quantities: form factors and values of the velocity transfer, 
uj = v ■ v' . Physical meson masses have been used as input where dimensional 
results are required. The only dependence on the lattice scale enters in expressing 
the heavy-light pseudoscalar masses in physical units for the 1/M extrapolations 
of the form factors and in the short distance logarithmic corrections of equa- 
tion (|ll]). The systematic error induced is much smaller than the other errors we 
will consider below. 

We believe our results show good evidence for light spectator mass indepen- 
dence of the form factors. We made this independence an assumption in order 
to use the results for the heaviest light spectator quark only: these results have 
smaller statistical error and are available for all our heavy quarks. In the previous 
section we argued that we could not rule out a 5% error in Ax (10% in channels 
involving momenta for the light meson) due to this assumption. The cor- 
responding errors in A2 and V were 20% and 10% respectively. In the region 
near q 2 ^^, dT/dq 2 is dominated by Ai and the contributions of V and A 2 are 
kinematically suppressed, being 5% or less for q 2 > 16 GeV 2 . The -> ^2 and 
1 — > y/2 channels explore lower q 2 values where the V and A 2 contributions are 
less suppressed, of order 15-20% (see figure H). If the spectator quark mass inde- 
pendence assumption were violated by as much as 20% for V and A 2 , the error in 
the differential decay rate would only be of the order of 2% for the non-\/2 chan- 
nels, and 6-8% for the v2 channels noted above, smaller in both cases than the 
likely 10% error induced by the same assumption for A\ (the differential decay 
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rate is proportional to the square of the form factors)^. Systematic errors due to 
this approximation are within our statistical errors. A conservative estimate of 
the systematic error in the differential decay rate arising from this source can be 
obtained by adding in quadrature the 10% error due to A\ and the worst-case 
8% error from V and A2, to get a final error of 13%. We will see below that 
discretisation errors are expected to be almost twice as large as this. 

We discussed in section ||] that we could not dismiss having 10% errors in the 
form factors arising from discretisation errors at the value of (3 used. In conse- 
quence, we will admit a 20% error in the decay rate which translates into a 10% 
error in \V u b\. These errors are the largest among the systematic effects consid- 
ered, and they almost entirely determine our overall systematic error. Combining 
errors from discretisation and our assumption of light quark spectator mass inde- 
pendence of the form factors in quadrature, we finally obtain an estimated 24% 
systematic error in the decay rate, which should be added to the results of table [| 
and equations (p6|)-(p8|). 



7 Conclusions 

We have presented a method for extracting the CKM matrix element \V u b\ from 
experimental measurements of the exclusive decay B° — > p + l~vi. We have shown 
how lattice simulations provide a model-independent framework in which nonper- 
turbative QCD corrections can be systematically incorporated in the evaluation 
of the differential decay rate dT(B° — > p + l~vi)/dq 2 in a region of q 2 values from 
the charm threshold up to the zero recoil point, q^ax- We have measured points 
spanning roughly two-thirds of this region of q 2 . One of the main features of 
the approach presented here is that it does not rely upon large, difficult and 
model-dependent extrapolations in q 2 from g^ ax down to values close to q 2 = 0, 
where the semileptonic decays are dominated by charmed final states, as in other 
previously suggested methods. 

Discounting experimental errors, the results of table [| will allow the determi- 
nation of the CKM matrix element \ V u b\ with a theoretical uncertainty of less than 
10% statistical and 12% systematic. Alternatively, determinations of \V u b\ (with 
similar theoretical ambiguities) can be obtained from the overall normalisation 
of the hadronic effects parametrised by a, found to be 

a = 4.6±g;| (stat) ± 0.6 (syst) GeV. (29) 

Both statistical and systematic theoretical errors will be reduced in forthcoming 
lattice simulations with smaller lattice spacing (to reduce discretisation errors, 
currently the principal systematic effect) and higher statistics. 



3 We have studied the effect on the results of the fit of equation (27) of having a 10% error 



in the v2 light meson momentum channels. This changes the parameter a by 5%. 
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Currently, only an upper bound for the total decay rate T(B° — > p + l~i>i) is 
known experimentally, but new results are forthcoming. It is important that the 
differential decay rate near g^ ax be measured, since, as we have shown, such a 
measurement can provide a clean determination of \V u b\ in the same way that 
similar measurements for the B — > D* semileptonic decay have successfully been 
used to extract an accurate value for |V^|. 

In this paper we have also determined the q 2 dependence of the form factor 
A\ (which dominates the decay rate close to the zero recoil point) for the B° — ► 
p + l~i>i decay in the region from q 2 = to q^ax- The form factors V and A 2 , for 
the same semileptonic decay, have been studied in a region of q 2 around g^ax- 

Finally, relations between radiative and semileptonic B decays into light vec- 
tor mesons, predicted by HQS, have been examined and found to hold within 
errors for the ratio A\/2iT 2 and within 20% for the ratio V/2T\ at the B scale. 

Using extra momentum channels where the heavy-light meson is not at rest, 
but u) is nevertheless nearly constant as the heavy quark mass changes []36|| , has 
been essential for studying q 2 , or equivalently 00, dependences of A x , Ai/2iT 2 , 
V/2Ti and dT/dq 2 . Knowledge of the q 2 dependence of dT/dq 2 has allowed us to 
determine the parameter a from a set of five measurements as well as from the 
zero recoil channel, giving us further confidence in our procedure for extracting 

|K&|- 
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